Optimal operation of a power plant

ABSTRACT

In a method and computer program product for optimal operation of a power plant and a power plant optimising system an optimisation module ( 1 ) minimises a cost function (J[u],J[P]) that comprises a deviation of lifetime of plant components (LTp(τ)) from a desired nominal lifetime trajectory (LTn(τ)). This is done by, at a given time, determining future values of input values (u(τ),P(τ)) such as control values (P(τ)) or process values (u(τ)) to the plant and simulating, in a simulation module ( 2 ), the behaviour of the plant up to a given future time. Corresponding lifetime values are determined in the simulation, and incorporated in an objective function. The optimisation module ( 1 ) minimises the cost function (J[u],J[P]) by varying the input values (u(τ),P(τ)). As a result, it is possible to operate the plant such that component lifetime (LTp(τ)) follows the desired trajectory (LTn(τ)).

FIELD OF THE INVENTION

[0001] The invention relates to the field of process control, in particular to a method and computer program for optimal operation of a power plant, and a power plant optimising system as described in the preamble of claims 1, 11 and 12, respectively.

BACKGROUND OF THE INVENTION

[0002] In known methods and systems for optimal operation of a power plant, plant operation is controlled by control values generated by a control system and is optimised in order to minimise fuel costs while achieving nominal required output values for produced power and process steam. This is done by, at a given time, determining future values of control values and simulating a behaviour of the plant up to a given future time. From the simulation, fuel costs are determined and an objective function comprising the fuel costs and costs for buying power from another source is computed. In order to determine optimal control values, the control values are varied and the simulation is repeated until a minimum of the objective function is found. The above procedure is also called “Unit commitment problem” (UCP). A mathematical treatment is given in “The Generalised Unit Commitment Problem”, Ross Baldick, IEEE Transactions on Power Systems, 10(1):465-475, February 1995. The optimisation is constrained by the fact that predetermined operating limits of the plant may not be exceeded.

[0003] “Optimal-Maintenance Modelling on Finite Time with Technology Replacement and Changing Repair Costs”, Jason W. Rupe, IEEE 2000 Proc. Annual Reliability and Maintainability Symposium, shows a method to determine when to replace or maintain components in a system with limited lifetime.

[0004] “Optimal Generation Scheduling with Ramping Costs”, C. Wang, and S. Shahidehpour, IEEE Trans. on Power Systems, Vol. 10, no. 1, pp 60-67. February 1995 , shows how to include an effect of fatigue stress and a resulting steam turbine rotor depreciation in scheduling algorithms for minimising operation cost.

DESCRIPTION OF THE INVENTION

[0005] It is an object of the invention to create a method and computer program for optimal operation of a power plant and a power plant optimising system of the type mentioned initially that allow a better combination of technical and economical considerations.

[0006] These objects are achieved by a method and computer program for optimal operation of a power plant and a power plant optimising system according to the claims 1, 11 and 12.

[0007] According to the invention, in a method and computer program for optimal operation of a power plant and a power plant optimising system according to the independent claims, a lifetime of plant components is determined when simulating a future behaviour of the plant, and a cost function to be minimised by the optimisation comprises a deviation of said lifetime from a nominal lifetime trajectory.

[0008] In this way, it is possible to prescribe a lifetime trajectory for a plurality of plant components and to operate the plant so that one or more components reach the end of their lifetime, or a state in which maintenance is necessary, at a predetermined time. It is therefore no longer necessary to schedule maintenance according to planned or expected plant operation, but plant operation can be modified in order to control when to perform maintenance.

[0009] In a preferred variant of the invention, plant process variables representing internal states of the plant are determined in the simulation and a difference between said plant process variables and nominal process variables is included in the cost function. The nominal process values are determined in order to reach given production and/or economic goals.

[0010] Further preferred embodiments are evident from the dependent patent claims.

BRIEF DESCRIPTION OF THE DRAWINGS

[0011] The subject matter of the invention will be explained in more detail in the following text with reference to preferred exemplary embodiments which are illustrated in the attached drawings, in which:

[0012]FIG. 1 shows a block diagram corresponding to the invention; and

[0013]FIG. 2 shows a block diagram corresponding to a preferred embodiment of the invention.

[0014] The reference symbols used in the drawings, and their meanings, are listed in summary form in the list of reference symbols. In principle, identical parts are provided with the same reference symbols in the figures.

DETAILED DESCRIPTION OF PREFERRED EMBODIMENTS

[0015]FIG. 1 shows a block diagram corresponding to the invention. An optimisation module 1 is provided with a nominal or desired lifetime trajectory LTn(τ) and with nominal or desired process variables Pn(τ). The optimisation module 1 generates power plant operating parameters or control values u(τ) that are input to a simulation module 2. The simulation module 2 comprises a simulation model of an actual power plant 3 and determines values for a power plant's internal states or process variables P(τ), and values of a lifetime trajectory LTp(τ) of plant components. The values determined by the simulation module 2 are input to the optimisation module 1. The values LTn(τ),LTp(τ),Pn(τ),P(τ),u(τ) mentioned above are, in general, vectors that change over time, i.e. each vector represents a trajectory in a multidimensional space. For example, LTp(τ) represents a lifetime trajectory for a plurality of pressure vessels, pipes, turbines, pumps etc. over the next hours, days or months. P(τ) represents e.g. steam pressures and temperatures and their gradients, or metal temperatures of critical components.

[0016] The control values u(τ) represent plant input or operation parameters that are used to control the plant, such as valve, pump or guide vane set points or actual positions, or set points for local controllers of pressure, level, temperature, fuel flow or steam mass flow. From a trajectory of control values u(τ) that is given for a time interval from a present time t to a future time t+T, where T is called “prediction horizon” or “optimisation horizon”, the simulation model determines a trajectory of the state of the power plant, i.e. the values of a plurality of variables that represent the state of the plant as it changes with time. From the state trajectory, trajectories LTp(τ),P(τ) for the values output by the simulation module 2 are determined for said time interval. In a preferred embodiment of the invention, the simulation model is a dynamic model of the power plant and its components.

[0017] The nominal process variables Pn(τ) are determined from process knowledge, taking into account a computer model of process physics and preferably optimising process efficiency and/or another optimisation criterion such as reliability. They are either predetermined or computed online during plant operation.

[0018] The lifetime trajectory LTp(τ) of plant components represents a remaining lifetime for a plurality of plant components. A lifetime value for a component may, depending on convention, either represent a “remaining lifetime” or a “consumed lifetime”. The latter are sometimes denoted as “equivalent operating hours” (EOH). Whichever convention is adopted does not affect the invention: A remaining lifetime decreases steadily from a total lifetime down to or beyond zero, while a consumed lifetime increases steadily up to and beyond the total lifetime. The total lifetime is the time span that the component is designed to be in operation under given nominal conditions. Lifetime boundary values determine when servicing or replacement of a component should take place—at the latest, when remaining lifetime is zero or less. After such servicing or replacement, the component's lifetime is reset accordingly.

[0019] Plant components are e.g. pressurised vessels, pipes or turbine blades. They are designed to withstand normal operating conditions for a given amount of time, i.e. for their total lifetime under nominal conditions. Subjecting a component to more severe operating conditions such as higher temperatures or temperature gradients increases component degradation and reduces component lifetime. This in turn increases long term costs for operating the power plant. Nevertheless it may be advantageous to determine nominal process values Pn(τ) that exceed normal operating limits and incur the corresponding operating losses if gains from production are sufficiently high. According to the invention, the corresponding change in lifetime is also taken into account and controlled by operating the plant accordingly.

[0020] The lifetime trajectory LTp(τ) is computed from models associated with various plant components, such as the following models, which are well know to one versed in the art:

[0021] Microcrack growth rate models for components operating at elevated temperatures, such as boiler pressure parts, steam pipes and steam turbine rotors. See for example “Life Assessment Technology for Fossil Power Plants”, R. Viswanathan, Proc. American Power Conf., Vol. 57-III: 1749-1768, 1995 and “A Lifetime Prediction Model for Thermomechanical Fatigue based on Microcrack Propagation”, M. P. Miller et.al, ASTM STP 1186, 1993, pp. 35-49.

[0022] Corrosion, oxidation, and erosion models based on extensions to microcrack growth models as described in “Corrosion Fatigue: Mechanics, Metallurgy, Electrochemistry, and Engineering”, T. Crooker and B. Leis editors, ASTM 1983.

[0023] Mechanical and thermomechanic stress models to determine stresses created in plant components, e.g. stresses in a steam turbine rotor shaft due to load oscillations and/or steam temperature changes.

[0024] A lifetime LT_(i)(t) is associated with an i'th component, for example, a pressurised vessel subject to constant and/or cyclic stresses. For such a vessel, LT_(i)(t) is, for example, computed from a crack growth model that determines an average current crack size a_(i)(t). The lifetime is determined as

LT_(i)(t)=k(a _(crit) −a _(i)(t))

[0025] where is a_(crit) a critical crack size and k is a normalising coefficient. When LT_(i)(t) reaches zero or becomes negative, the component is assumed to have reached or exceeded its lifetime.

[0026] Any other method for calculating component lifetime from operating conditions may be used, such as the ones shown in the following publications of Deutscher Dampfkessel-Ausschuβ, “TRD 300: Calculation of boiler strength”, “TRD 301: Cylindrical shells under internal pressure”, “TRD 508 Annex 1: Additional tests on components—methods for the calculation of components having time dependent design strength values”, 1978-79.

[0027] In the inventive method, the following steps are performed

[0028] a) The optimisation module 1 determines the trajectory of plant control values u(τ) from the present time t for a predetermined duration T and transmits said plant control values u(τ)to the simulation module 2. The first time that this step is executed, plant control values are e.g. set to constant values for the duration T or are determined from stored standard trajectories corresponding to a given operating condition such as a plant start-up or a load change. Later executions of this step preferably use a variation of the plant control values u(τ) of an earlier step.

[0029] b) The plant model incorporated in the simulation module 2 determines from the trajectory of control values u(τ), by simulation of the power plant behaviour, a trajectory of the state of the power plant. From said state, corresponding trajectories of process variables P(τ) are determined. Control values are, for example, valve positions, motor speeds, electric power flow, etc. Process variables P(τ) represent the physical process of main interest. In the case of a power plant this is the water/steam cycle with process variables such as mass flows, temperatures, pressures, fluid levels etc. Some process variables P(τ) are identical to states of the power plant in the control sense, others are computed from the states by linear or non-linear and/or dynamic submodels. From the process variables P(τ), simulated values of lifetime trajectory LTp(τ) of plant components are computed.

[0030] c) A cost function J[u] is computed as an integral of an objective function that comprises a deviation of the lifetime trajectory LTp(τ) from the given desired lifetime trajectory LTn(τ). In a preferred embodiment of the invention, this cost function J[u] is defined as the functional ${{J\lbrack u\rbrack} = {\int\limits_{t}^{t + T}\left. ||{{R_{LT}(\tau)}\left( {{{LTn}(\tau)} - {{LTp}(\tau)}} \right)}||{\tau} \right.}},$

[0031] where R_(LT) is a weighting matrix. It determines a relative influence of the magnitude of components ΔLT_(i)(t) of the difference between nominal and simulated actual lifetime vectors LTn(τ),LTp(τ). The norm ∥ . . . ∥ used in the integral is any norm, for example a 2-norm or an infinity norm.

[0032] In a preferred embodiment of the invention, the above objective function also includes a deviation of the process variables P(τ) from the desired nominal process variables Pn(τ). In this embodiment, the cost function J[u] is preferably defined as ${{J\lbrack u\rbrack} = {\int\limits_{t}^{t + T}\left. ||{{R_{LT}(\tau)}\left( {{{LTn}(\tau)} - {{LTp}(\tau)}} \right)}||{+ \left. ||{{R_{p}(\tau)}\left( {{{Pn}(\tau)} - {P(\tau)}} \right)}||{\tau} \right.} \right.}},$

[0033] where R_(p) is a corresponding weighting matrix for process variables. It allows prescription of the relative importance of deviations of individual variables with respect to one another and to lifetime values. Both R_(LT) and R_(p) are time dependent and positive semidefinite. If a 2-norm is used, the function is ${J\lbrack u\rbrack} = {\int\limits_{t}^{t + T}\left. ||{{R_{LT}(\tau)}\left( {{{LTn}(\tau)} - {{LTp}(\tau)}} \right)}||{}_{2}{+ \left. ||{{R_{p}(\tau)}\left( {{{Pn}(\tau)} - {P(\tau)}} \right)}||{}_{2}{{\tau}.} \right.} \right.}$

[0034] d) Steps a) through c) are repeated iteratively, with the optimisation module varying the trajectory of plant control values u: [t,t+T]→R^(n) until an optimised lower value of the cost function J[u] is arrived at As the trajectories LTp(τ),P(τ) depend upon plant control values u(τ), the optimisation routine finds a trajectory of control values u(τ) or strategy u*: [t,t+T]→R^(n) that minimises the cost function J[u]. Appropriate optimisation and simulation methods are well known from the field of model predictive control (MPC).

[0035] In another preferred embodiment of the invention, the cost function in the steps shown above is expressed as J[P] in terms of the process variables, e.g. as ${J\lbrack P\rbrack} = {\int\limits_{t}^{t + T}\left. ||{{R_{LT}(\tau)}\left( {{{LTn}(\tau)} - {{LTp}(\tau)}} \right)}||{+ \left. ||{{R_{p}(\tau)}\left( {{{Pn}(\tau)} - {P(\tau)}} \right)}||{\tau} \right.} \right.}$

[0036] and the optimisation according to step d) determines

P*=argminJ[P]

[0037] with the strategy P*: [t,t+T]→R^(q) minimising the cost function J[P]. In this embodiment, the optimisation determines the optimal trajectory P*(τ) of process variables. The task of determining corresponding control values u(τ) that achieve the optimal trajectory P*(τ) of process variables is left achieved by a control function or control system for of the plant. This function or system may be considered to reside, in the structure according to the figures, in either the simulation (2) or in the optimisation (1) module.

[0038] In the case in which the cost function J[u] is minimised the plant control values u(τ) are considered to be input values that are varied for the simulation and by optimisation. When the cost function J[P] is minimised, then the process variables P(τ) are considered to be said input values.

[0039] In a further preferred embodiment of the invention, the cost function incorporates a term that represents a difference between the lifetime vector LT_(p)(T_(E)) and the nominal lifetime vector LTn(T_(E)) at a given predetermined end time T_(E), weighted with an appropriate weighting matrix R_(E). For example, the cost function J[u] or J[P] is $\int\limits_{t}^{t + T}\left. ||{{R_{LT}(\tau)}\left( {{{LTn}(\tau)} - {{LTp}(\tau)}} \right)}||{+ \left. ||{{R_{p}(\tau)}\left( {{{Pn}(\tau)} - {P(\tau)}} \right)}||{{\tau} +}||{R_{E}\left( {{{LTn}\left( T_{E} \right)} - {{LTp}\left( T_{E} \right)}} \right)}||. \right.} \right.$

[0040] This term realises a soft constraint that penalises not reaching the nominal lifetime vector LTn(T_(E)) at time T_(E).

[0041] In another preferred embodiment of the invention, this penalty is achieved by the time dependency of R_(LT), i.e. by having relatively large entries in R_(LT) at the time T_(E). If following the exact trajectory of the nominal lifetime vector LTn(τ) is not important most of the time, the entries in R_(LT) are set to relatively small values for these times.

[0042] In yet a further preferred embodiment of the invention, the optimisation is subject to hard constraints on one or more of the plant control values u(τ) and/or of the process variables P(τ).

[0043] In another preferred embodiment of the invention, the nominal lifetime vector LTn(τ) is determined such that the nominal lifetime consumption of components is co-ordinated, with a plurality of given components reaching an associated predetermined lifetime limit value at the same time. For example, if maintenance of a given component such as a boiler involves shutting down other components such as pumps, tanks etc. simultaneously, then maintenance of these other components may take place when the boiler is shut down. Accordingly, the nominal lifetime of these other components is scheduled to reach or at least allowed by choice of weighting values to reach an associated lifetime limit value at essentially the same time as the boiler.

[0044] The features of the above preferred embodiments can be varied and combined in a manner obvious to one versed in the art. For example, one embodiment comprises a cost function J[P] using infinity-norms and incorporating a soft constraint on the lifetime vector LTp(T_(E)) at the end time T_(E).

[0045] The optimisation step c) involves the problem of how to find the optimal solution efficiently, since the search space determined by the plurality of plant control values and the duration T is large. Standard optimisation procedures based on e.g. dynamic programming or the Hamilton-Jacobi method may be used, or dedicated solutions as described, for example, in “The Generalised Unit Commitment Problem”, Ross Baldick, IEEE Transactions on Power Systems, 10(1):465-475, February 1995. If the cost function incorporates a hard constraint on the lifetime at the end time, then a solution to the corresponding optimisation problem is found e.g. by the Pontryaguin maximum principle. Alternative optimisation routines based on brute force methods coupled with complexity reducing heuristics are also possible, for example, simulated annealing and genetic algorithms.

[0046] In a preferred embodiment of the invention, at time t, the control strategy, i.e. one or more of the future values u: [t,t+T]→R^(n) corresponding to optimal plant operation and/or corresponding values of the optimal trajectory P*(τ) of process variables are displayed graphically to an operator, providing her with guidance for controlling the plant.

[0047]FIG. 2 shows a block diagram corresponding to another preferred embodiment of the invention, in which, at time t, an optimal control value u*(t) is applied to the actual power plant 3, resulting in automatic optimised plant operation. Preferably, measured plant values y(t) of the plant are input to the simulation module 2 in order to update the simulation model. Appropriate modelling and model update techniques are well known from the fields of optimal control and model predictive control. In particular, receding horizon control techniques give useful results. The power plant is, for example, a gas turbine, a combined cycle power plant, or of a boiler-turbine-generator type.

[0048] A computer program product according to the invention comprises a computer readable medium, having thereon: computer program code means to make, when said program is loaded in a computer that is operationally connected to a power plant, the computer execute the method according to the invention. A computer program according to the invention comprises said computer program code.

[0049] A power plant optimising system for optimal operation of a power plant according to the invention comprises data processing means such as control computers and hardware and/or software modules, including:

[0050] a) a simulation module 2 that is configured to determine, from a trajectory of input values, by simulation of the power plant's behaviour, corresponding trajectories of process variables P(τ), and

[0051] b) an optimisation module 1 that is configured to determine a plurality of trajectories of input values and to determine a value of a cost function as an integral of an objective function, which objective function comprises at least one trajectory determined by the simulation module 2.

[0052] The simulation module 2 is configured to determine, from the trajectory of input values, a trajectory representing a lifetime LTp(τ) of plant components, and the objective function comprises a difference of said lifetime LTp(τ) from a nominal lifetime trajectory LTn(τ).

[0053] In a preferred embodiment of the power plant optimising system, the objective function comprises a difference between the trajectory of process variables P(τ) and a trajectory of associated nominal process variables Pn(τ). A further preferred embodiment of the power plant optimising system comprises means for displaying at least part of the trajectory of input values to an operator and/or means for using at least part of said values to control the actual plant (3).

[0054] The invention allows the operator or an automatic controller to incorporate not only short-term financial considerations into a plant operating strategy, but also allows scheduling and consideration of the influence of lifetime consumption on components, in order to meet and/or optimise planned replacement and maintenance schedules.

List of Designations

[0055]1 optimisation module

[0056]2 simulation module

[0057]3 actual power plant

[0058] u(τ) control values

[0059] u*(t) optimal control values

[0060] LTp(τ) lifetime trajectory

[0061] LTn(τ) desired or nominal lifetime trajectory

[0062] P(τ) process variables

[0063] Pn(τ) desired or nominal process variables

[0064] y(t) measured plant values

[0065] R_(1LT),R_(P),R_(E) weight matrices 

1. A method for optimal operation of a power plant, where the method comprises the steps of a) determining a trajectory of input values (u(τ),P(τ)) b) determining, from the trajectory of input values (u(τ),P(τ)), corresponding trajectories of process variables P(τ), c) computing a cost function (J[u],J[P]) as an integral of an objective function that comprises at least one trajectory determined in the simulation, d) iteratively repeating steps a) through c) with an optimisation module (1) varying the trajectory of input values (u(τ),P(τ)) until an optimised lower value of the cost function (J[u],J[P]) is arrived at, characterized in that in step b), a trajectory representing a lifetime LTp(τ) of plant components is determined from the trajectory of input values (u(τ),P(τ)), and that, in step c), the objective function comprises a difference of said lifetime LTp(τ) from a nominal lifetime trajectory LTn(τ).
 2. Method according to claim 1, wherein the input values are control values u(τ), the trajectories of process variables P(τ) are determined from u(τ) by simulation of the power plant's behaviour, and the optimisation step finds a trajectory of control values u(τ) that minimises the cost function J[u].
 3. Method according to claim 1, wherein the input values are the process variables P(τ) and the optimisation step finds a trajectory of process variables P(τ) that minimises the cost function J[P].
 4. Method according to one of the preceding claims, wherein in step c) the objective function comprises a difference between the trajectory of process variables P(τ) and a trajectory of associated nominal process variables Pn(τ).
 5. Method according to one of the preceding claims, wherein in step c) the cost function incorporates a term that represents a difference between the lifetime vector LTp(T_(E)) and the nominal lifetime vector LTn(T_(E)) at a given predetermined end time T_(E).
 6. Method according to one of the preceding claims, wherein in step c) the optimisation is subject to hard constraints on one or more of the plant control values u(τ) and/or of the process variables P(τ).
 7. Method according to one of the preceding claims, characterized in that at least part of the trajectory of input values (u(τ),P(τ)) corresponding to the optimised cost is displayed to an operator.
 8. Method according to one of the preceding claims, characterized in that at least part of the trajectory of input values (u(τ),P(τ)) corresponding to the optimised cost is used to control the actual plant (3).
 9. Method according to one of the preceding claims, wherein the cost function (J[u],J[P]) comprises the term ${\int\limits_{t}^{t + T}\left. ||{{R_{LT}(\tau)}\left( {{{LTn}(\tau)} - {{LTp}(\tau)}} \right)}||{\tau} \right.},$

where R_(LT) is a weighting matrix, t is a present time and T is a predetermined duration.
 10. Method according to claim 9, wherein the cost function (J[u],J[P]) comprises the term ${\int\limits_{t}^{t + T}\left. ||{{R_{LT}(\tau)}\left( {{{LTn}(\tau)} - {{LTp}(\tau)}} \right)}||{+ \left. ||{{R_{p}(\tau)}\left( {{{Pn}(\tau)} - {P(\tau)}} \right)}||{\tau} \right.} \right.},$

where R_(1LT),R_(P) are weighting matrices, t is a present time and T is a predetermined duration.
 11. A computer program comprising computer program code means to make, when said program is loaded in a computer that is operationally connected to a power plant, the computer execute the method according to one of claims 1 to
 10. 12. A power plant optimising system for optimal operation of a power plant that comprises a) a simulation module (2) that is configured to determine, from a trajectory of input values (u(τ),P(τ)), by simulation of the power plant's behaviour, corresponding trajectories of process variables P(τ), b) an optimisation module (1) that is configured to determine a plurality of trajectories of input values (u(τ),P(τ)) and to determine a value of a cost function (J[u],J[P]) as an integral of an objective function, which objective function comprises at least one trajectory determined by the simulation module (2), characterized in that the simulation module (2) is configured to determine, from the trajectory of input values (u(τ),P(τ)), a trajectory representing a lifetime LTp(τ) of plant components, and that the objective function comprises a difference of said lifetime LTp(τ) from a nominal lifetime trajectory LTn(τ). 